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Turbulent penetration can occur when hot and cold fluids mix in a horizontal T-junction pipe at nuclear plants. 
Caused by the unstable turbulent penetration, temperature fluctuations with large amplitude and high frequency 
can lead to time-varying wall thermal stress and even thermal fatigue on the inner wall. Numerous cases, however, 
exist where inner wall temperatures cannot be measured and only outer wall temperature measurements are feasi- 
ble. Therefore, it is one of the popular research areas in nuclear science and engineering to estimate temperature 
fluctuations on the inner wall from measurements of outer wall temperatures without damaging the structure of 
the pipe. In this study, both the one-dimensional (1D) and the two-dimensional (2D) inverse heat conduction 
problem (IHCP) were solved to estimate the temperature fluctuations on the inner wall. First, numerical models of 
both the 1D and the 2D direct heat conduction problem (DHCP) were structured in MATLAB, based on the finite 
difference method with an implicit scheme. Second, both the 1D IHCP and the 2D IHCP were solved by the 
steepest descent method (SDM), and the DHCP results of temperatures on the outer wall were used to estimate 
the temperature fluctuations on the inner wall. Third, we compared the temperature fluctuations on the inner wall 
estimated by the 1D IHCP with those estimated by the 2D IHCP in four cases: (1) when the maximum disturbance 
of temperature of fluid inside the pipe was 3°C, (2) when the maximum disturbance of temperature of fluid inside 


the pipe was 30°C, (3) when the maximum disturbance of temperature of fluid inside the pipe was 160°C, and (4) 


when the fluid temperatures inside the pipe were random from 50°C to 210°C. 
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Introduction 


Thermal fatigue is one of the most important degrada- 
tion issues for heavy nuclear components. Fatigue dam- 
age evaluation by measuring the outer wall temperatures 
is useful for maintaining the integrity of heavy nuclear 
components. Inverse heat conduction problem (IHCP) is 
often applied in those engineering problems in which 
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direct measurements of temperatures on the inner wall 
are not possible. 

The solution of IHCP has received considerable inter- 
est because of numerous important applications in 
science and engineering. In a recent paper Dou et al. [1] 
investigated an IHCP to determine the heat transfer cha- 
racteristics of a water jet on a stainless steel plate. Liu 
and Takase [2] solved an IHCP to obtain the surface heat 
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flux and surface temperature of a block by using the 
measured inner temperature data as input. 

The solution of IHCP is not new, particularly in the 
field to estimate the thermal fatigue of a pipe under 
thermal stratification. Lu et al. [3, 4] investigated both 
two-dimensional and three-dimensional IHCPs to esti- 
mate the unknown transient fluid temperatures near the 
inner wall of a pipe elbow. However, the solution of 
IHCP is still new, in the field to estimate the thermal fa- 
tigue of a pipe under turbulent penetration. 

A variety of numerical techniques have been proposed 
for the solution of IHCPs: (1) Levenberg-Marquardt me- 
thod (LMM) (An IHCP was solved by LMM to identify 
the thermophysical properties of soil [5]. An IHCP was 
solved by LMM to determine the optimal geometry of 
the fillers between two conductive bodies in a two- 
dimensional multiple-region domain [6]). (2) Least- 
squares method (LSM) (Using the weighted LSM, we 
formulated an inverse problem of textile thickness-heat 
conductivity-porosity determination (IPT (THP) D) into a 
function minimization problem [7]. The temperature- 
dependent thermal diffusivities of two alloys were deter- 
mined by solving an IHCP by LSM [8]). (3) Conjugate 
gradient method (CGM) (An IHCP was solved by CGM 
to identify the heat losses at the jaws of a tensile testing 
machine [9]. An IHCP was solved by CGM to estimate 
the time-dependent heat flux by using temperature dis- 
tribution at a point in a three-layer system [10]). (4) Ge- 
netic algorithm (GA) (An IHCP was solved by GA to 
identify the conductive and radiative parameters of a se- 
mitransparent sample [11]. An IHCP was solved by GA 
to determine the thermal conductivity of lumber [12]). (5) 
Steepest descent method (SDM) (An IHCP was solved by 
the SDM to identify the boundary of unknown inclusions 
inside an object by applying a heat source and measuring 
the induced temperature near the boundary of the sample 
[13]. A solution of an IHCP was obtained by the SDM to 
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determine the waste heat flux from a helicon plasma dis- 
charge by using the transient surface temperature mea- 
surements obtained from infrared thermography [14]). 
The SDM, which has great potential in solving IHCP, 
transforms the inverse problem to the solution of three 
problems, namely, the direct problem, the sensitivity 
problem, and the adjoint problem [15]. Although the 
SDM has great potential in solving IHCP, its validity in 
solving IHCP for a T-junction pipe under turbulent pene- 
tration has never been examined. 

In this study, the SDM was used to solve both the one- 
dimensional (1D) and the two-dimensional (2D) IHCPs 
in order to estimate the inner wall temperature fluctua- 
tions in a horizontal T-junction pipe under turbulent pe- 
netration. 


Physical and mathematical models 


(1) Direct heat conduction problem 

The fatigue test specimen is horizontal and shown in 
Fig. 1. Turbulent penetration can occur when hot and 
cold fluids mix in a horizontal T-junction pipe. Our work 
is aimed to estimate fatigue of pipes under turbulent pe- 
netration. Fig. 2 shows the two-dimensional physical 
model of the main pipe and the branch pipe of a horizon- 
tal T-junction pipe. The inner and outer walls are both 
subjected to the convective boundary condition. 

The governing equation for the direct heat conduction 
problem (DHCP) to calculate the temperature fluctua- 
tions on the inner wall is 

2 2 
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where T is the temperature; æ is the thermal diffusivity 
of the stainless steel pipe; x and y are the space variables; 
T is the time variable; and © is the computational do- 
main. 


T upstream 


Fig. 1 Fatigue test specimen. 
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Fig. 2 Physical model. 


The corresponding boundary conditions are 
1D oo = houl ik E Tanl (x,y) € outer wall (2) 
Aua = hin (Toin 5 Ton) (x,y) einner wall (3) 
where / is the thermal conductivity of the stainless 
steel pipe, n is the direction vector normal to the boun- 
dary wall, his the convection heat transfer coefficient, Te 
is the fluid temperature, and 7,, is the wall temperature. 
The subscripts in and out represent the inner and 
outer walls, respectively. 
The initial condition is 
T(x, y,0) = 1) (x,y) (4) 
Where 7 (x, y) is the initial pipe temperature. 
(2) Inverse heat conduction problem 
When we know all the boundary conditions and the 
initial condition, the unsteady two-dimensional DHCP 
can be solved with commercial CFD codes such as 
FLUENT. However, numerous cases exist where the 
boundary condition on the inner wall is unknown. In 
these cases, a direct solution with commercial CFD code 
is not available, and we have to develop programs to 
solve an IHCP 
During the IHCP solution, a DHCP numerical model 
is structured with an implicit scheme based on the finite 
difference method. Then the boundary condition on the 
inner wall is repeatedly assumed and the corresponding 
DHCP is repeatedly solved to compare the calculated 
outer wall temperatures with the measured values. The 
IHCP can be treated as an optimization problem; there- 
fore, we use an objective function based on the least- 
squares method to solve the IHCP: 


N M 
J= > » [Tamea ~ ere li (5) 


n=] m=1 
where Tn, m, cai 18 the calculated temperature on the out- 
er wall, Ta, m, mea is the measured temperature on the outer 
wall, m denotes the number of measure points on the 
outer wall, n denotes the number of time steps, M is the 
total number of measure points on the outer wall, and N 
is the total number of time steps. 
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The consistency between the calculated temperatures 
and the measured temperatures acts as a constraint on the 
optimization problem. During the optimization, the 
steepest descent method (SDM) is used to determine the 
search direction and step size for each iteration. The 
temperature on the outer wall varies with the temperature 
on the inner wall; i.e., Tp, m, cai is a function of Tp: 

Irme = oln) (6) 

where 7), m, cai is the calculated temperature on the out- 

er wall, and Tin is the temperature on the inner wall. Fur- 

thermore, J is a function of Tin, and we can get the opti- 
mization problem: 

minJ =J (Tp) (7) 

The optimization problem can be solved by the fol- 
lowing iterative formula: 

Ta? = Fy + Aya ®© 
where b denotes the number of iterations, Ta® is the 
approximate root of the optimization problem after b 
iterations, d” is the search direction when Tn equals 7; AOA 
and ^; is the step size when T;n equals Te ) When we use 

the SDM, d® will be the steepest descent direction: 
d® =-VI (Ta) (9) 

and /, will meet the following condition: 
I(T + Ad) = min J (Ti? + Ad) (10) 
>0 


The computational procedure for the solution of the 
IHCP is summarized as follows and shown in Fig. 3: 


Initial 7, 


eo 
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M 
break | J 
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Fig. 3 Flowchart of solution of IHCP. 


Step1. Set an initial 7; pD, use a small positive number 
gas the stopping criterion, and set b equal to 1; 
Step2. Solve the DHCP; 


GUO Zhouchao et al. 


Step3. Compute the steepest descent direction d®; 


Step4. If d” meets the condition: la a < €, terminate 


the iteration process. Otherwise, compute the step size 4%; 
Step5. Compute p {ED set b equal to b+1, and go to 
Step2. 


Results and discussion 


Both the 1D IHCP and the 2D IHCP were solved with 
temperatures on the outer wall obtained from the DHCP 
results to estimate the temperatures on the inner wall. 
Furthermore, we compared the temperatures on the inner 
wall estimated by 1D IHCP with those estimated by 2D 
IHCP in four cases: (1) when the maximum disturbance 
of temperature of fluid inside the pipe was 3°C, (2) when 
the maximum disturbance of temperature of fluid inside 
the pipe was 30°C, (3) when the maximum disturbance of 
temperature of fluid inside the pipe was 160°C, and (4) 
when the fluid temperatures inside the pipe were random 
from 50°C to 210°C. 

(1) When the maximum disturbance of temperature of 
fluid inside the pipe was 3 °C 

Fig. 4 shows the comparisons at 0°, 90°, 180°, and 
270° from 0 s to 15 s, and we can see that the 1D IHCP 
and 2D IHCP results show the same trend. Table 1 shows 
the average, maximum, and minimum relative errors, and 
we can see that the maximum relative error of 0.32% 
(shown in red and bold) occurs at 180°. 
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Fig. 4 Comparisons between 1D and 2D estimated inner wall 
temperatures under case (1)(when the maximum dis- 
turbance of temperature of fluid inside the pipe was 
3°C). 


Table 1 Average, maximum, and minimum relative errors 
under case (1). 


Ave (%) Max (%) Min (%) 
IHCP_1D (0°) 0.08 0.20 0.00 
IHCP_1D (90°) 0.13 0.31 0.01 
IHCP_1D (180°) 0.11 0.32 0.00 
IHCP_1D (270°) 0.19 0.31 0.02 
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(2) When the maximum disturbance of temperature of 
fluid inside the pipe was 30°C 

Fig. 5 shows the comparisons at 0°, 90°, 180°, and 
270° from 0 s to 15 s, and we can see that the 1D IHCP 
and 2D IHCP results also show the same trend. The av- 
erage, maximum, and minimum relative errors are shown 
in Table 2. The maximum relative error of 1.74% (shown 
in red and bold) also occurs at 180°. But the average, 
maximum, and minimum relative errors increase as the 
maximum disturbance of temperature of fluid inside the 
pipe increases from 3°C to 30°C. 


Tin / °C 


Tin / °C 


Fig.5 Comparisons between 1D and 2D estimated inner wall 
temperatures under case (2)(when the maximum dis- 
turbance of temperature of fluid inside the pipe was 


30°C). 


Table 2 Average, maximum, and minimum relative errors 
under case (2). 


Ave (%) Max (%) Min (%) 
THCP_1D (0°) 1.02 1.45 0.68 
THCP_1D (90°) 0.50 1.42 0.07 
THCP_1D (180°) 1.04 1.74 0.36 
THCP_1D (270°) 0.26 0.95 0.01 


(3) When the maximum disturbance of temperature of 
fluid inside the pipe was 160°C 

Fig. 6 shows the comparisons at 0°, 90°, 180°, and 
270° from 0 s to 15 s, and we can see that the 1D IHCP 
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and 2D IHCP results also show the same trend. The av- 
erage, maximum, and minimum relative errors are shown 
in Table 3. The maximum relative error of 7.29% (shown 
in red and bold) also occurs at 180°. But the average, 
maximum, and minimum relative errors increase as the 
maximum disturbance of temperature of fluid inside the 
pipe increases from 30°C to 160°C. 
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Fig. 6 Comparisons between 1D and 2D estimated inner wall 
temperatures under case (3)(when the maximum dis- 
turbance of temperature of fluid inside the pipe was 
160°C). 


Table 3 Average, maximum, and minimum relative errors 
under case (3). 


Ave (%) Max (%) Min (%) 
IHCP_1D (0°) 3.42 5.15 1.93 
IHCP_1D (90°) 2.30 5.91 0.29 
THCP_1D (180°) 4.68 7.29 1.83 
IHCP_1D (270°) 1.43 4.97 0.08 


(4) When the fluid temperatures inside the pipe were 
random from 50°C to 210°C 

Fig.7 shows the comparisons at 0°, 90°, 180°, and 
270° from 0 s to 15 s, and we can see that the 1D IHCP 
and 2D IHCP results also show the same trend. The av- 
erage, maximum, and minimum relative errors are shown 
in Table 4. The maximum relative error of 7.87% (shown 
in red and bold) occurs at 90°. But the average, maxi- 
mum, and minimum relative errors under case (4) are 
larger than those under case (3). 
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Fig. 7 Comparisons between 1D and 2D estimated inner wall 
temperatures under case (4)(when the fluid tempera- 
tures inside the pipe were random from 50°C to 
210°C). 


Table 4 Average, maximum, and minimum relative errors 
under case (4). 


Ave (%) Max (%) Min (%) 
THCP_1D (0°) 2.93 5.29 1.29 
THCP_1D (90°) 4.80 7.87 2.03 
THCP_1D (180°) 2.23 3.44 0.39 
THCP_1D (270°) 2.12 4.01 0.12 


Conclusions 


Both the 1D IHCP program and the 2D IHCP program 
were developed in MATLAB to estimate the inner wall 
temperature fluctuations in a horizontal T-junction pipe 
under turbulent penetration. Furthermore, we compared 
the 1D IHCP results with the 2D IHCP results under 
some specific boundary conditions on the inner wall. The 
results and comparisons showed that 

(1) an algorithm based on the SDM can be applied to 
solve the IHCP to estimate the inner wall temperature 
fluctuations in a horizontal T-junction pipe under turbu- 
lent penetration; 

(2) the average relative errors in all four cases above 
were less than 5%, indicating that under some specific 
cases the 1D IHCP can be used to estimate the inner wall 
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temperature fluctuations for online monitoring of thermal 
fatigue in a mixing T-junction pipe. 
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